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Abstract. The effects of general relativity (GR) on the hydrodynamics and 
neutrino transport are examined during the critical shock reheating phase 
of core collapse supernovae. We find that core collapse computed with GR 
hydrodynamics results in a substantially more compact core structure out 
to the shock, the shock radius at stagnation being reduced by a factor of 2. 
The inflow speed of material behind the shock is also increased by a factor 
of 2 throughout most of the evolution. We have developed a code for general 
relativistic multigroup flux-limited diffusion (MGFLD) in static spacetimes 
and compared the steady-state neutrino distributions for selected time slices 
of post-bounce models with those computed with Newtonian MGFLD. The 
GR transport calculations show the expected reductions in neutrino lumi¬ 
nosities and rms energies from redshift and curvature effects. Although the 
effects of GR on the hydrodynamics and neutrino transport seem to work 
against shock revival, the core configurations are sufficiently different that 
no firm conclusions can be drawn, except that simulations of core collapse 
supernovae using Newtonian hydrodynamics and transport are not realistic. 


1. Introduction 

General relativity (GR) is an essential component in the realistic modeling of core collapse 
supernovae because of the very strong gravitational fields in the vicinity of the collapsed 
core of a star. Hydrodynamics and neutrino transport are closely connected in this 
problem, and as we will show, GR can have a profound effect on each of these, especially 
in the critical phase of shock reheating. The detection of neutrinos from supernova 1987A 
(Bionta et al. 1987, Hirata et al. 1987) and the hope of detecting neutrino signatures 
from future supernovae, with next-generation detectors, is additional motivation for an 
accurate general relativistic treatment of the neutrino transport in numerical simulations. 
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acknowledge the hospitality and support of the Institute for Theoretical Physics in Santa Barbara, which 
is supported by the National Science Foundation under contract PHY94-07194. 
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Radii vs. Post-Bounce Time 

S15s7b 



Post-Bounce Time, sec 

Figure 1. Shock and gain radii vs. post-bounce time for model S15s7b. Both cases are 
calculated with Newtonian radiation transport. 


We have developed a code for general relativistic multigroup flux-limited diffusion 
(MGFLD) that computes the neutrino transport in a static background metric. This 
is the hrst step in the development of a fully GR MGFLD code. The GR metric used 
is ds"^ = a^c^dt^ — h'^dr'^ — r‘^{d9‘^ + Odcp"^). This metric is allowed to evolve in the 
hydrodynamics calculations, but is “frozen” in the transport calculations, which are 
then performed treating the metric as static. Two precollapse models, a ISM© model 
(S15s7b) and a 25Mq model (S25s7b) (Woosley & Weaver 1995; Weaver & Woosley 
1997), were evolved through core collapse, bounce, and to approximately 800 ms after 
bounce in three sets of simulations. The hrst was with Newtonian hydrodynamics and 
Newtonian radiation transport, the second was with GR hydrodynamics and Newtonian 
transport, and the third was with both GR hydrodynamics and transport. In addition to 
these three sets of simulations, stationary-state neutrino distributions were computed for 
various post-bounce time slices in these models in order to isolate the effects associated 
with each of the metric components. 
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RMS Energy vs. Radius 


15 Solar Mass Model (S15s7b) 



Figure 2. RMS energy vs. radius for the in model S15s7b. This hgure contrasts 
stationary state GR and Newtonian transport at tpb = 114 ms. 

2. The Role of General Relativity 

The effects of GR are seen quite clearly in the hydrodynamic evolution of the initial 
models. GR hydrodynamics produces a much more compact post-bounce structure than 
Newtonian hydrodynamics. After fpb = 0.4 s in both models (fpb being the post-bounce 
time), the radius of the shock and the gain radii are reduced by a factor of 2 in the GR 
calculations, as shown in Figure 1. 

Also strongly affected by GR is the flow velocity between the shock and the proto¬ 
neutron star. Because matter falls through a greater potential well to reach the shock 
in the GR calculation, GR preshock and therefore postshock velocities are larger than 
their Newtonian counterparts, again by a factor of approximately 2. 

The main effect of GR on the neutrino rms energies is the redshift of the neutrinos 
after they decouple from the matter, which is governed by the metric parameter a. The 
rms energies are reduced by a factor of a evaluated at the //-sphere. A smaller effect is 
a slight outward shift of the //-sphere resulting from the (non-unity) value of the metric 
parameter b~^, which causes the neutrinos to decouple outside the //-sphere at a lower 
temperature. These effects are shown in Figure 2. Also shown are the independent 
effects of a and b on the neutrino transport. 
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Luminosity vs. Radius 

15 Solar Mass Model (S15s7b) 



Figure 3. Luminosity vs. radius for the VeS in model S15s7b. This hgure contrasts 
stationary state GR and Newtonian transport at tpb = 114 ms. 

GR reduces the neutrino luminosity by three effects: redshift, governed by the 
metric parameter a; the difference in the local clock rates at the emission surface and 
the observer radius, also governed by the metric parameter a; and the reduction of the 
neutrino flux, governed by the metric parameter b~^. All three of these effects are of 
roughly equal magnitude and reduce the luminosities by a total factor of a?b~^ evaluated 
at the //-sphere, as shown in Figure 3. This hgure also shows the independent effects of 
a and b on the stationary-state neutrino transport. 

The shock heating rate is proportional to the product of the luminosity and the 
square of the z/g rms energy. This means that any percentage change in these quantities 
add together. Looking at Figures 2 and 3, we see a decrease of 8% in the Ue luminosity 
and a 3% reduction in the z/g rms energy. These combine to give a 14% reduction in the 
heating rate. Similar results are obtained for the z/g’s. 

These are signihcant differences [e.g., see Burrows & Goshy (1993), Janka & Muller 
(1996), Mezzacappa et al. (1998), and Messer et al. (1998)] and serve to illustrate the 
point that modeling core collapse supernovae without GR hydrodynamics and transport 
leads to results that cannot be interpreted as realistic. For more information on this 
subject, the reader is referred to Bruenn, De Nisco, & Mezzacappa (1998). 
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RMS Energy vs. Post-Bounce Time 
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Figure 4. RMS energy vs. post-bounce time for z/e’s and z/e’s in model S15s7b. 

3. Nucleosynthesis 

r-process nucleosynthesis is believed to occur in a neutrino-driven wind emanating from 
the proto-neutron star after the successful launch of the shock. The r-process yields are 
a function of the z/e and z/e luminosities and rms energies. The luminosities affect the 
entropy, mass loss rate, and expansion time scale associated with the wind, and the rms 
energies determine the neutronization of the wind. As an example of the impact of GR 
on the r-process, consider its effect on the rms energies. Because the z/g-sphere lies below 
the zze-sphere, the z/e’s suffer a greater emergent redshift. As can be seen in Figure 4, 
general relativistic transport and hydrodynamics affects the ratio of the z/e and z/e rms 
energies. This differential redshifting affects, in turn, the ratio of the number of z/g’s 
to z/e’s, at a given neutrino energy, and therefore, the neutronization of the wind. This 
suggests that general relativistic hydrodynamics and transport will be required to obtain 
accurate r-process yields. 
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